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Sulphur abundances in halo stars from Multiplet 3 at 1045 nm* 
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Sulphur is a volatile a-element which is not locked into dust grains in the interstellar medium (ISM). Hence, its abundance 
does not need to be corrected for dust depletion when comparing the ISM to the stellar atmospheres. The abundance of 
sulphur in the photosphere of metal-poor stars is a matter of debate: according to some authors, [S/Fe] versus [Fe/H] forms 
a plateau at low metallicity, while, according to other studies, there is a large scatter or perhaps a bimodal distribution. In 
metal-poor stars sulphur is detectable by its lines of Mult. 1 at 920 nm, but this range is heavily contaminated by telluric 
absorptions, and one line of the multiplet is blended by the hydrogen Paschen £ line. We study the possibility of using 
Mult. 3 (at 1045 nm) for deriving the sulphur abundance because this range, now observable at the VLT with the infra-red 
spectrograph CRIRES, is little contaminated by telluric absorption and not affected by blends at least in metal-poor stars. 
We compare the abundances derived from Multiplets 1 and 3, taking into account NLTE corrections and 3D effects. Here 
we present the results for a sample of four stars, although the scatter is less pronounced than in previous analysis, we 
cannot find a plateau in [S/Fe], and confirm the scatter of the sulphur abundance at low metallicity. 
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I Introduction 

The light elements between O and Ti of even atomic number 
are referred to as a-elements because they are mainly pro- 
duced by successively adding an a-particle from nucleus to 
nucleus. Among them is sulphur, a volatile element which 
is not locked into dust in the interstellar medium (ISM), so 
that the sulphur abundance derived for the ISM can be di- 
rectly compared to the sulphur abundance derived in stars. 

The a-elements are crucial probes of the chemical evo- 
lution of a stellar population: they are almost exclusively re- 
leased by Type II Supernovae, while the iron peak elements 
are produced by Type II SNe, but also, in large amounts, 
by Type la Supernovae. Progenitors of Type II and Type la 
supernovae have very different lifetimes, making the abun- 
dance ratio of a-elements to iron-peak elements a powerful 
diagnostic of the chemical evolution and star formation his- 
tory of a galaxy. In the Milky Way, stars of lower metallicity 
are characterised by higher a to iron abundance ratios than 
found in the Sun and stars of solar metallicity. This is usu- 
ally interpreted in terms of the lower contribution of Type 
la SNe. Systems which are characterised by low or bursting 
star formation, like dwarf Spheroidal galaxies, give time to 
Type la SNe to explode before the enrichment due to Type 

II SNe has greatly increased. Consequently such systems 
display rather low a to iron ratios even at low metallicities. 



Thus the abundance of a-elements is an important property 
of any stellar population. For the study of chemical evolu- 
tion in external galaxies, the more readily available objects 
are Blue Compact galaxies (BCGs) through analysis of the 
emission line spectra, and Damped Ly— a systems (DLAs) 
through the analysis of resonance absorption lines. In both 
groups of objects, sulphur is relatively easy to measure in 
the form of ISM emission (for the BCGs) or absorption (for 
the DLAs) lines. 

The investigation of sulphur abundances in the stel- 
lar photospheres started with the pioneering work of 



Wall erstein & Contil (119641) . They determined the sulphur 



abundance in six stars out of a sample of nine stars in 
six Ga lactic Clusters. Later on. IClegg. Tomkin. & Lambert 
dl98ll) determined the su lphur content in 20 F- G -type stars 
with [Fe/H]> -1, while [Francois! dl987l [l988h studied 13 
and 12 metal -poor stars, respectively. 

Not many sulphur lines are available in the observed 
stellar spectra. There is a forbidden line from the ground 
level of S I at 1082.1 nm, which is weak, blended, but for 
which lower and upper level populations are very close 
to local thermal equilibrium (LTE). This absorption line i s 



measurable in the spectrum of solar-like stars (Ryde , 2006). 



* Using data from CRIRES at the ESO-VLT, Programme 079.D-0434. 
** Corresponding author: e-mail: Elisabetta.Caffau@obspm.fr 



However, it becomes undetectable below [Fe/H]< —0.5. 
The permitted Mult. 8 at 675 nm is weak, but not blended, 
neither contaminated by telluric absorption. For these tran- 
sitions the assumption that high and low levels are in LTE 
is a good approxima tion, and the same holds for the lines 
of Mult. 6 at 869 nm dTakada-Hidai et alll2002h . Both Mul- 
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tiplets 6 and 8 are weak, as a consequence detectable only 
in stars of solar, or moderately sub-solar metallicity (down 
to [Fe/H]~ —1.5). Below such metallicity, usually Mult. 1 
at 920 nm is used, which is, on the other hand, contami- 
nated by telluric absorption. Moreover one of the lines of 
the triplet (922.8 nm) is located in the blue wing of the hy- 
drogen Paschen £ line, and the presence of the hydrogen line 
makes the abundance determination more difficult. 

The infra-red lines of Mult. 3 at 1045 nm are well suited 
to measure the sulphur abundance. These lines are not as 
strong as the components of Mult. 1, but easily detectable in 
very metal-poor stars. No blend has to be taken into account 
in the case of metal-poor stars because the only extra line 
present in the range, the Fe I line at 1045.5 nm blending the 
strongest S I line of the triplet, vanishes at sub-solar metal- 
licity. Telluric absorptions in this range are less abundant 
and less strong than in the case of Mult. 1 . Observing the 
1045 nm sulphur lines provides thus a possibility to obtain 
a reliable sulphur abundance in very metal-poor stars. 

Caffau et al. (2005) analysed the sulphur abundance in 



Table 2 S I lines of Mult. 3 



a sample of Galactic stars. They investigated, when avail- 
able, the lines of Multiplets 6, 8, and 1. Because of the 
blend of the 922.8 nm S I line of Mult. 1 with the hydro- 
gen Paschen £, they performed line profile fitting to derive 
the sulphur content for all the lines available. The line pro- 
file fitting procedure permitted to reproduce the line pro- 
file of the hydrogen. In their study they have suggested that 
in the range -2.5 <[Fe/H]< -2.0 the [S/Fe] ratio shows 
either a large scatter or a bimodal behaviour. Most of the 
stars lie on a "plateau" at about [S/Fe]=+0.4, while a non 
negligible number of stars shows a "high" value of [S/Fe], 
around +0.8. This behaviour has no proposed interpreta- 
tion, but it may be supposed to be due to systematic errors 
affecti ng only the analys is of Mult. 1 . In fact, in the sam- 
ple of ICaffau et al.l (12005b . the determination of [S/Fe] in 
this range of metallicity is based mainly on the non con- 
taminated lines of Mult. 1 . On the other hand, Nissen et al. 



(2004J) and iNissen et al.l (120071) derived the sulphur abun- 
dance from equivalent width (EW) measurements of the 
lines of multiplets 8 and 1, and they find a plateau at low 
metallicity in the [S/Fe] versus [Fe/H] plot, with no sign of 
bimodal distribution or scatter. It is therefore of great inter- 
est to verify this puzzling finding by the use of an indepen- 
dent and, probably better, diagnostic of the sulphur abun- 
dances, as can be afforded by the 1045 nm lines. 

2 Model atmospheres and atomic data 

This analysis was performed by means of ID hydrostatic 
model atm ospheres. We t hen a pplied NLTE corrections 
taken from Takeda et al.l (2005) and 3D corrections de- 
rived by means of CO 5 BO LP 3D hydrodynamical model s 



dFrevtag. Steffen. & Dorchi 120021; IWedemever et aUl2003l) . 
The ID model atmospheres used to derive the sulphur abun- 
dance were obtained from the Linu x version of ATLAS, as 
described in ISbordone et al.l (12.004ft . To derive the 1D-LTE 



A 

(nm) air 


Transition 


log gf 


Xlo 

(eV) 


1045.5449 
1045.6757 
1045.9406 


3go _ 3p 2 
3 S? - 3 P 

3go _ 3pi 


0.26 
-0.43 
0.04 


6.86 
6.86 
6.86 



abundance of sulphur we measured the equivalent widths 
(EW) with the IRAF task splot and derived abundances 



Kurucz, 



through WI DTH (for details see ICastelliL 12005 
Il993ll2005l) . 

The 3D-C0 5 BOLD models used to compute the 3D- 
corre ctions (listed in Tab le Q]) are taken from the CIFIST 
grid dLudwig et al. , 2009). For each 3D model, we used 
as reference model the ID model obtained as the hor- 
izontal average of each 3D snapshot over surfaces of 
equal (Rosseland) optical depth, and the ID hydrostatic 
model computed with the LHD code with the same atmo- 
spheric parameter as the 3D model. Both these ID models 
share the micro-physics and radiative transfer scheme with 
CQ 5 BOL D. A de s criptio n of these models can be found in 



Caffau & Ludwia (120071). We comput ed 3D c orrections as 



describ ed in ICaffau & Ludwrgl (2007) and in Caffau et al. 
(1201 Oh 



The 3D corrections are given in Table [4] in the case 
of a micro-turbulence of 1.5km s _1 used in the ID ref- 
erence model. AW(£ m i c ) = A(S) 3D — A(S)^ 3D ^ in col- 
umn 7 takes into account the effects of granulation, mean- 
ing the influence that the horizontal fluctuations around the 
mean stratification have on the abundance determination. 
A (2) (£mic,aMLT) =A(S) 3D -A(S) 



;iD, 



in column 8 ac- 
counts for the total 3D effect, meaning both the effect of 
horizontal fluctuations and the different mean temperature 
structure due to the different treatment of convection in 3D 
hydrodynamical models and in ID mixing-length models. 
The atomic data of the S I lines of Mult. 3, the same as 

inlj 



Caffau et al.l (I2007h . are summarised in Tabled 



3 Analysis of CRIRES data 

We observed four metal-poor dwarf stars with CRIRES 
dKaeufl et all 12004). The observations, programme 079. D- 
0434 (E. Caffau), have been carried out in service mode, 
for details see Table [3] The spectra were observed using 
detector integration times of 30 s or 10 s (for the brighter 
stars), the total integration time for each star is given in Ta- 
ble[3] The observation template employed, required nodding 
along the slit, with an amplitude of 10". The AO star HIP 
28910 was also observed, to check the position of telluric 
lines. 

The analysis was conducted using the extracted, com- 
bined and calibrated spectra provided by ESO. We com- 
pared this data with that obtained for one star by running 
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Table 1 Hydrodynamical model atmospheres 



Model 


Teff 


logg 


[Fe/H] 


At 






Box-size 
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h 






km 3 






d3t50g35mmlOn01 


4930 


3.5 


-1.0 


117 


59732.2 


X 


59732.2 


X 


30199.5 


d3t50g35mm20n01 


4980 


3.5 


-2.0 


151 


59732.3 


X 


59732.3 


X 


30199.5 


d3t50g45mml0n03 


5060 


4.5 


-1.0 


18 


5084.46 


X 


5084.46 


X 


2477.64 


d3t59g35mm20n01 


5860 


3.5 


-2.0 


31 


59308.6 


X 


59308.6 


X 


38658.3 


d3t59g35mm30n01 


5870 


3.5 


-3.0 


31 


59308.6 


X 


59308.6 


X 


38658.3 


d3t59g45mm20n01 


5920 


4.5 


-2.0 


7 


6020.0 


X 


6020.0 


X 


3841.66 



Table 3 S I lines of Mult. 3 



columns 2-4 of Table @] In column 5 the 1D-LTE abun- 
dance with the unce rtainty related to th e EW measurement 



Star 


J 


date 


total integration time 




mag 




s 


BD-05°3640 


7.62 


24-07-2007 


2x 2160 


HD 140283 


6.01 


24-07-2007 


300 


HD 165195 


4.89 


24-07-2007 


400 


HD 181743 


8.62 


19-06-2007 


2160 






23-07-2007 


2160 


HD 21 1998 


4.15 


22-07-2007 


160 



is given. While in ICaffau et al. (120051) no NLTE correc 



tion was appli ed, here NLTE corr ections are interpolated 
in the table of iTakeda et al. (2005) and are listed in col- 
umn 6 of Table |4] The 3D corrections for each star are 
given in column 7 and 8. In column 9 [S/Fe] is given; 
both NLTE and A' 2 ) 3D corrections are applied. We adopt 
A(S)=7.16 for the s olar sulphur abundance, as described in 
Caffau et alJ (bold). In column 10 the sulphur abundance 



a set of reduction routines based on the ESO CRIRES 
pipeline 1.6.0, but did not find any significant difference 
with respect to the reduced data provided by ESO. 

Our sample consisted of four Galactic stars from 



from lCaffau et alJ d2005 ) is reported, a nd in column 1 1 ther e 
is the relative NLTE correction from ITakeda et alJ ([2005). 
The comparison in [S/Fe] is not straightforward. To derive 
the [S/Fe] given in column 9 of Table |4j both NLTE and 
3D-A 1 - 2 ' corrections are applied. As one can see in the Ta- 
ble, these two corrections roughly cancel each other, so that 
the IDlte abundance we derive is close to the abundance 
obtained after the application of NL TE and 3D correctio ns. 
Another difference with respect to ICaffau etaL d2005l) is 
the solar reference abundance , A(S)q = 7.16 in this analy- 
dance derived is the one in Tab leffl H D 181743 has also sis. whilelCaffau et all J2005I) adopted A(S) = 7.21 from 



Caffau et al 



the sample of 
BD-05°3640, HD 181743 



(120051) . Three of them, 
and HD 211998, have been 
analysed by ICaffau et alJ (12005b . and the sulphur abun- 



been analysed by iNissen et al.1 (|2004). As stellar param- lAnders & Grevessd (1989). For this reasons we suggest to 



eter T eff /logg [Fe/H] they determined 5863 K74.32/-1.93, 
and their sulphur abundance was A(S)=5.61, in agree- 
ment, within erro rs with our result.. HP 21 1998 has also 
been analysed by Israelian & Rebolol ( 2001 ) and Francois 



compare the A(S) 1D-LTE obtained in this work and in 



Caffau et alJ (12005b . in column 4 and 10 of Table |4]respec- 



dl987l) . Ilsraelian & Rebolol feOOllT derived A(S)=5.80 with 



the stellar parameters T eff /log g/[Fe/H] of 5271K/3.36/- 
1.25. The difference of more than 0.3 dex with respect to 
our result is anyway compatible within errors, in fact for this 
star their error is about 0.3 dex. We stress also a difference 
of 0.3 dex in met allicity. Ou r analy sis is in close agreement 
with the one of iFrancoisI Jl987l) . who find s A(S)=6.20 
For HP 140283 the values in Tableware from lNissen et al.1 
d2004l) . 

Our goal was to understand if Mult. 3 can be used in or- 
der to derive the sulphur abundance in metal-poor stars. It 
could be th at the large scatter in [S/Fe] at low metallicity 
described in Caffau et al. (2005) is due to systematic effects 
arising from the use of one specific S I multiplet for abun- 
dance determination. A sample of 4 stars is insufficient to 
firmly establish whether this is the case or not, but it can 
provide useful insights to direct further research. 

The stel lar parameter s for th e stars of our sample are 
taken from ICaffau et alJ (120051) and are summarised in 



tively. 

For st ar BP — 05°3640, t he "high" sulphur abundance 
derived in lCaffau et al.1 ([2005) is here confirmed. The differ- 
ence of 0. 16 dex between the two analysis is large but within 
the errors. The same is true for HP 21 1998: for this star 
the NLTE correction for Mult. 1 is 0.07 dex larger than for 
Mult. 3. For HP 140283 the two analyses are in good agree- 
ment. In Fig.Q]we show the observed spectrum (solid-black 
line) of HP 140283 compared to an ATLAS+SYNTHE syn- 
thetic spectrum (dashed-red line). For star HP 181743 we 
find a difference b etween this analysis and the results of 



Caffau etal J (12005) larger than expected from the error bars. 



The previous analysis relied on Mult. 6, the observed spec- 
trum was observed with UVES@VLT, with resolution of 
43 000, dichroic #2 and cross disperser #4. In this spectrum 
a line was clearly visible at the wavelength of the Mult. 6, 
with a FWHM compatible with the iron line close by. The 



spectr um showed no fringing as suggested by INissen et al. 



(2007). A possible explanation for the presence of this line 
is that it is a contribution of the sky, not perfectly subtracted. 
In fact at the time of that observation the Moon was bright 
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Table 4 Stellar parameters and sulphur abundances. To derive [S/Fe] in column 9, NLTE and 3D corrections have been 
applied. For the solar reference we adopt A(S) Q = 7.16. 



Star 


Teff 


logg 


[Fe/H] 


A(S)-1DLTE 


Alte 


A w 


AW 


[S/Fe] 


A(s) uos 


Alte 




K 








Mult. 3 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


BD-05°3640 


5020 


4.61 


-1.19 


6.61 ±0.10 


-0.02 


0.03 


0.05 


0.67 ±0.14 


6.77 ±0.10 


-0.03 


HD 140283 


5690 


3.69 


-2.42 


5.12 ±0.05 


-0.15 


0.00 


0.08 


0.31 ±0.11 


5.10 ±0.10 


-0.16 


HD 181743 


5970 


4.40 


-1.81 


5.73 ±0.10 


-0.07 


0.02 


0.13 


0.44 ±0.14 


5.83 u ±0.10 


-0.04 


HD 211998 


5210 


3.36 


-1.56 


6.15 ±0.10 


-0.12 


0.06 


0.12 


0.55 ±0.14 


6.33 ±0.10 


-0.19 



AW(e mic ) = A(S) 3D - J 4(S) (3D) 

A< 2 > (Uc «mlt) = A(S) 3T) - A(S) 
u: upp er limit from this w ork 
C05: lcaffauetai]fc005h 
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Fig. 1 The spectral region of the S I lines of Mult. 3 in 
the case of the metal-poor star HD 140283 (solid black) 
over-imposed to an ATLAS+SYNTHE synthetic spectrum 
(5690K/3.69/2.5, A(S)=5.1 1, dashed red). 



(68% illumination) and close to the star (29° distance). Due 
to the large discrepancy we found with the analysis of the 
lines of Mult. 3, we reduced the archival UVES spectra of 
this star taken during ESO programmes 67.B-0474(A) (P. I. 
Holmberg) as well as 67.D-0106(A) (P.I. Nissen). In the first 
programme a single 600', s, R~60000 spectrum was pro- 
duced with CD #3, while in the second case a set of 400 s, 
R~60000 spectra were produced with image slicer #1 and 
CD #2. With these data we cannot detect the Mult. 6 and 
are only able to give an upper limit for the sulphur abun- 
dance, which is reported i n Table [4] We confirm the re- 



sult of Nissen et al 



d2007l) about the Mult. 6, and find a 
sulphur abundance from Mult. 3 very close to their value. 
In fact, according to their LTE analysis, [S/Fe]=0.29 with 
A(S) Q = 7.20, while we find [S/Fe]= 0.38 ± 0.010 with 
A(S) = 7.16. 



4 Discussion 



Caffau et alj d2005l) derived the sulphur abundance in a sam- 
ple of metal-poor stars and included in their analysis the re- 
sults available in the literature for a large number of Galactic 
stars. Their complete sample spans a wide range in metal- 
licities, from solar to 1/1000 of the solar metallicity. As vis- 



ible in their figure 10, the S abundance, as a function of the 
stellar iron content, has a negative trend for stars of solar 
metallicity or for slightly metal-poor stars (—1.0 < [Fe/H] < 
±0.5). For low metallicity ([Fe/H]< —1.5) a large scat- 
ter or a bimodal behaviour is visible. This trend, which at 
the moment has no physical explanation, could have been 
attributed either to the fact that the sulphur abundance in 
metal-poor stars was mainly derived from the S I lines of 
Mult. 1 where telluric absorption can be a problem, or to the 
neglect of NLTE effects. NLTE effects are not very large, 
about -0. 1 dex, and in part they are compensated by 3D ef- 
fects. The scatter we find in this work from Mult. 3 from 
our 4 stars is larger than 0.3 dex, compatible only within 2a 
with a plateau. In Fig. HI our 1D-L TE results are compared 
to the results o f lCaffau et al. (2005), which have been renor- 
malised to A(S) Q = 7.16. 

Caffau etail d2007l) investigated the sulphur abundance 
derived from the S I lines of Mult. 3 in the Sun and a small 
sample of solar metallicity stars. They compared this abun- 
dance with the one obtained from the weaker lines of Multi- 
plets 6 and 8. Their analysis was carried out before CRIRES 
was available and the needed spectral range was obtained 
by forcing the extraction of the last order in the UVES 
860 nm standard setting of a few bright, solar-metallicity 
stars taken from the Paranal Observatory Project dataset 



(Bag nulo et all 120031) . Except for the Sun, the abundances 



obtained from the different lines of S I are in agreement. 



Taka da-Hidai et al.l d2002h derived sulphur abundances 
in six giants and 61 dwarfs, in the range — 3 < [Fe/H] < 
±0.5, using Mult. 6, and found a linear increase of [S/Fe] 
versus [Fe/H] with decreasing metallicity. 



Taka da-Hidai et al.l (120051) analysed Multiplets 6 and 1 



in 21 metal-poor stars. They find a good agreement of A(S) 
from the two multiplets for [Fe/H]> —2, while for the most 
metal-poor stars the abundances from Mult. 6 are larger than 
from Mult. 1. In their analysis [S/Fe] as a function of the 
metallicity forms a plateau. 

Nissen et all d2007b observed Mult. 3 of G 29-23, us- 



ing CRIRES science verification data, and obtained a good 
agreement of the LTE sulphur abundance with the re- 
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Fig. 2 The 1D-LTE results of ICaffau et all (120051) . from 
their analysis (asterisk s) and taken from the literature 
(cross es ), meaning f rom Chen et al. (2002); Ecuvillon et al 
d2004l): llsraelian & Rebolol d200lh: iNissen et al.1 d2C 
iRvde & Lambertl d2004 :~lTakada-Hid ai et all foom . 
compared to our value (filled circles) from Mult. 3. 
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Fig. 3 The 1D-LTE results of INissen et al.1 ( 120071) 
(crosses) are compared to our value (filled circles) from 
Mult. 3. 



the previous analysis relied on Mult. 1. We have to notice 
that for three stars the sulphur abundance we obtain from 
the lines of Mult. 3 is lower with respect to the previous 
analysis based on Mult. 1. For HD 181743, whose previous 
analysis was based on Mult. 6, the difference of 0.5 dex is 
above the expected error, and can now be explained with 
a contamina t ion fro m the sky in the spectrum analysed by 
ICaffauet al.ld2005b . 

From the results of this small sample we cannot confirm 
the existence of a plateau in the [S/H],[Fe/H] plane. The 
scatter is larger with respect to the estimated uncertainties, 
but a similar sc atter is presen t in [ Mg/Fe] versus [Fe/H] in 
the analysis of IGratton et al.l (120031) . The systematic error 
related to the oscillator strength of the lines of Mult. 1 can 
be neglected because it acts in the same way for the com- 
plete sample of stars. One has to take into account the un- 
certainties related to the temperature, but these error cannot 
remove the great scatter in [S/Fe]. In fact for the cooler star 
a change of 100K in the effective temperature is translated 
in a change of 0.10 dex in the sulphur abundance, while for 
the hotter star the same difference in temperature would pro- 
duce only 0.04 dex change in A(S). 

An analysis on Mult. 3 in an extended sample of stars 
would be useful to investigate the behaviour of [S/Fe] at 
low metallicity. 
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suit obtained from the other lines. In Fig. [3] the sample of 
Nisse n et al. (2007), scaled by 0.04 dex for the difference in 



the solar sulphur between the two analysis, is compared to 
our results. The two most metal-ri ch stars we st u died show 
a [S/Fe] larger than the values in INissen et al. (2007) for 
similar metallicity, while the other two stars fall on the av- 
erage. A scatter of about 0.3 dex in [S/Fe] for metal-poor 
stars is comparable to the scatter in [Mg/Fe] available in 
IGratton et al.l (12003b for a similar range in metallicity and 
in lAndrievskv et al.l (120101) for more metal-poor stars. 



5 Conclusions 



We did an analysis on the S I lines of Mult. 3 for four m etal - 
poor dwarfs stars of the sample of ICaffau et al.l d2005). The 
results of three out of the four candidates are consistent with 
the previous analysis within errors, and for these three stars 
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